Experimental determination of the quasi-projectile mass with measured neutrons 
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Abstract 



"The investigation of the isospin dependence of multifragmentation reactions relies on precise reconstruction of the 
fragmenting source. The criteria used to assign free emitted neutrons, detected with the TAMU Neutron Ball, to the 
^ 'quasi-projectile source are investigated in the framework of two different simulation codes. Overall and source-specific 
i O .detection efficiencies for multifragmentation events are found to be model independent. The equivalence of the two 
different methods used to assign experimentally detected charged particles and neutrons to the emitting source is shown. 
.The method used experimentally to determine quasi-projectile emitted free neutron multiplicity is found to be reasonably 
'accurate and sufficiently precise as to allow for the study of well-defined quasi-projectile sources. 
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1. Introduction 

The study of multifragmentation as a function of isospin 
has recently become an area of interest |ll-fl3jj because of 
its connection with the isospin-dependent part of the equa- 
tion of state 14{. Much of the early experimental work 
was performed by varying the isospin content, N/Z, of 
the beam and/or target and measuring obscrvablcs as a 
function of the N/Z of the reacting system. Later, several 
works showed that the key parameter in understanding 
isospin dynamics might be the isospin content of the frag- 
menting system rather than that of the reacting system 
@, El; 0, EH ■ For instance, the excited quasi-projectile 
(QP, the momentarly existing primary fragment, remnant 
of the projectile after a non-central collision) produced in 
peripheral and semi-peripheral collisions at intermediate 
energies can access a rather broad distribution in N/Z and 
the centroid of the distribution is shifted toward the valley 
of stability 0, [H, EH • More specifically, the importance 
of charge and mass characterization of the frag menting 
source was pointed out in several works Il6l | investi- 
gating the symmetry energy term of the nuclear equation 
of state. Isoscaling [lo[ and m-scaling [l7| have been shown 
to significantly improve when selecting quasi-projectiles 
within a narrow window of N/Z (l3l . [l6j|. Similarly, an 
agreement among different methods of extracting the sym- 
metry energy coefficient was reported in fl3j , provided that 
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stringent constraints were placed on the isotopic asymme- 
try of the QP source. 

The abundant emission of charged particles and neu- 
trons accompanying these reactions often tend to obscure 
the primary processes, making it difficult to characterize 
intermediate systems formed in the collision and study 
the underlying reaction mechanisms. Different approaches 
were used to determine the mean N/Z of the QP source. 
One approach uses isotopically resolved charged fragments 
to reconstruct the QP on an event-by-event basis and re- 
lies on models to include the evaporated neutrons in the 
measured N/Z of the QP [HOI, Il9j |. A second approach, 
which relies on the hypothesis of early fragment formation 
*20l . 21 1, uses fragment yield ratios of pairs of isobars to 



determine the average N/Z of the QP [H, S 

The availability of neutron calorimeters coupled to An 
charged particle detectors allows for direct event-by-event 
measurement of free emitted neutrons and charged parti- 
cles. This makes it possible to perform an event-by-event 
QP charge and mass measurements, independent of mod- 
els. 

In this work we will present and explore the reliability 
of the criteria that we use to experimentally assign de- 
tected free neutrons to the QP fragmenting source. For 
this purpose, overall and source-specific neutron detection 
efficiencies of the Neutron Ball detector will be calculated 
using the HIPSE-SIMON simulation code for reac- 

tion systems corresponding to the experimental data de- 
scribed in Sec. Neutrons will be assigned to the QP 
using both a source tag provided by the HIPSE code and 
a velocity selection experimentally used for charged parti- 
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cles (see Sec. 12. 2[) . The results obtained from both meth- 
ods of QP neutron assignment and their comparison will 
be reported throughout this paper to show their equiva- 
lence. In addition results obtained with the Constrained 
Molecular Dynamics model [25| for the same reaction sys- 
tems will be presented alongside HIPSE-SIMON results to 
illustrate the model independence of our criteria. Finally, 
the free neutron assignment method will be applied to our 
experimental data and the QP neutron distributions will 
be shown for three different systems. 

This article is organized as follows. In Sec. [2]we present 
the experimental setup and the QP reconstruction proce- 
dure using information on the emitted free neutrons. In 
Sec. |31 the calculation of the overall, QP and QT effi- 
ciencies for multifragmentation events is presented, as ob- 
tained by analysing two different models. Moreover, the 
reliability of the QP reconstruction as well as the effect of 
the experimental geometrical and detection efficiencies of 
the Neutron Ball are analysed. An application of the pre- 
sented criteria to experimental data is presented in Sec. 0] 
Conclusions are presented in Sec. [31 

2. Experimental conditions and procedures 

Several experiments were recently carried out at the 
Texas A&M University on the K500 superconducting Cy- 
clotron to investigate the isospin-dependent part of the 
equation of state by looking at QP fragmentation events. 
Such studies rely on the availability of an experimental 
setup which provides cvent-by-event information on the 
emitted charged particles and neutrons, and on a QP re- 
construction procedure capable of identifying the fragment- 
ing QP in charge, mass and exitation energy. In the fol- 
lowing we will describe the experimental setup and the QP 
reconstruction procedure used in these experiments. 

2.1. The experimental setup 

The 4tt NIMROD-ISiS array [H, ^ is used for the 
detection of charged particles and neutrons. The detector 
telescopes, arranged on 15 rings centered on the beam axis, 
are composed of one or two silicon detectors backed by a 
CsI(Tl) crystal with PMT readout. Isotopic resolution, 
essential for isospin physics analyses, is obtained up to 
Z = 17. Details on the mass identification resolution and 
procedures can be found in [28|, [2§| . The charged-particle 
array is housed inside the TAMU Neutron Ball (2?| , which 
measures the free neutron multiplicity. 

The Neutron Ball consists of 6 scintillator tanks, filled 
with a pseudocumene based liquid scintillator (NE224) 
loaded with 0.3 wt. % Gd which surrounds the vacuum 
chamber. Four large quadrant-shaped detectors are sand- 
wiched between the two end cups tanks, creating a cylin- 
der. Each scintillator tank is read by 3 or 4 fast photo- 
multiplier tubes (quadrants and cmisphcrcs, respectively). 
Following a nuclear reaction, prompt 7 rays and neutrons 
dissipate energy in the scintillaor, producing a flash. Rel- 
atively low energy neutrons tend to be thermalized in the 



scintillator and diffuse through the detector for a period of 
microseconds before undergoing capture reactions mainly 
with Gd and H nuclei. After the fast flash PMT signals, 
a delayed 100 pis gate is opened. The number of delayed 
flashes in the gate is proportional to the number of neu- 
trons emitted in the triggering reaction. Following this 
gate, a second 100 ps gate is opened. The multiplicity 
measured in the second gate provides information on the 
background. The beam is prevented from entering the de- 
tector during the counting gates by a deflector. The back- 
ground measurement also takes into account the possible 
activation of material in the chamber. 

Further details on the Neutron Ball geometry and elec- 
tronics, the detection mechanisms and the detector intrin- 
sic efficiency can be found in 26l. 30l]. 



2. 2. Quasi-projectile reconstruction procedure 

Particle and event selections are performed to select 
quasi-projectile (QP) fragmentation events. 

The quasi-projectile source is reconstructed on an event- 
by-event basis by associating charged particle fragments 
with the QP source based on a velocity selection in which 
each particle's longitudinal velocity (v z ) is compared to 
the v z of the largest fragment. For Z = 1 fragments, v z 
must be within ±65% of the largest fragment's v z . For 
Z = 2 and Z > 3, the cut is at ±60% and ±45%, respec- 
tively. This cut, later referred to as Vcut, is intended to 
remove fragments from non-projectile-likc sources and can 
be applied only to measured charged particles since the 
Neutron Ball does not measure neutron velocities. 

In order to select on QP fragmentation events, the to- 
tal charge of the fragments detected in the 4-7T array and 
included in the reconstruction was chosen to be between 
83% and 100% of the beam charge (SumZ). Finally, limits 
are placed on the deformation of the source, as measured 
by the quadrupole momentum, to select a class of events 
that are, on average, spherical [l2l ]. The quadrupole mo- 
mentum, calculated from the measured particle momenta 
in the quasi-projectile frame, X^Pf -/ Y^ii \p\ ■■> ^ s required 
to be between 0.4 and 2.7. This will be later referred to 
as Qcut. 

Free neutrons measured by the Neutron Ball were used 
to calculate the number of free neutrons emitted by the 
quasi-projectile (Nqp) using the relation [l2l ]: 

N QP = ^| ■ (1) 



n7 £ QT 



where Ndet is the measured neutron multiplicity. In order 
to calculate the number of free neutrons emitted by the 
QP, the total multiplicity is corrected by the efficiency for 
detecting free neutrons emitted from the quasi-projectile 
{sqp) and the quasi-target(egT) for our reactions. Equa- 
tion [T] assumes that QP and QT have the same relative 
N/Z as the projectile and target and that the main con- 
tributions to neutron emission come from QP and QT 
emissions, while other sources (such as, for instance, pre- 
equilibrium emission) are neglected. Np and Nt are the 
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projectile and target neutron numbers, respectively. The 
determination of the efficiencies for the QP and QT sources 
will be presented in Sec. 13.11 Through this formulation, 
the number of neutrons emitted from the QP can be cal- 
culated from only a total multiplicity of free neutrons. 

3. Analysis 

Recently, experimental data on 64 Zn+ 64 Zn, 70 Zn+ 70 Zn 
and 64 Ni+ 64 Ni at 35AMeV were used to investigate the 
isospin dependence of the caloric curve [3l[ and to compare 
different methods to extract infomation on the isospin- 
dependent part of the equation of state [l3| . Both studies 
pointed out the need for placing stringent constraints on 
the isotopic asymmetry of the analysed QP source. The 
uncertainty in the QP identification arises mainly from the 
assignment of detected free neutrons to the QP source. In 
order to develop a procedure for assigning free emitted 
neutrons to the QP source, a systematic study was per- 
formed on the predictions of two different models, HIPSE- 
SIMON [HI Hif and CoMD [25|, for the 64 Zn+ 64 Zn, 70 Zn 
and 64 Ni+ 64 Ni reactions. The number of events gener- 
ated with the HIPSE-SIMON code and the CoMD code is 
1600000 and 4040JL respectively. 

The HIPSE [23[ model is a phenomenological event 
generator dedicated to the description of nuclear collisions 
in the intermediate energy range. The approach is mainly 
based on the geometrical participant-spectator picture and 
accounts for both statistical and dynamical effects. The 
after-burn er p hase is achieved using the SIMON event 
generator 24| which takes into account all possible decay 
channels from neutron evaporation to symmetric fission. 
Information on the source of the final particle is provided 
in the output of the calculation which makes the code very 
suitable to investigate the source-specific efficiency of the 
Neutron Ball. 

The CoMD model [25| was also used, and the results 
were compared to the HIPSE-SIMON results. CoMD is 
a constrained molecular dynamics model for a fcrmionic 
system, suitable for the study of heavy colliding systems. 
The time evolution of each particle was followed, in our 
calculation, up to 3000 fm/c, to assure a full de-excitation 
of the fragment. No information on the particle-emitting 
source is provided by the code. 

An experimental filter was implemented to simulate the 
Neutron Ball response. The Neutron Ball efficiency was 
calculated in Ref.flU using the GCALOR code coupled 
to the GEANT-3 simulation package to simulate the Neu- 
tron Ball incorporated into NIMROD. Indeed, the large 
amount of material inside the Neutron Ball chamber, due 
to the presence of the charged particle array, can result 
in the scattering, absorption and generation of neutrons. 
The calculations were done for neutrons emitted isotrop- 
ically from the target with energies from 5 to 200 MeV 
(in 5 MeV steps). An efficiency of 0.6 was obtained for 
all neutron energies averaged over all angles. However, in 
heavy ion reactions around Fermi energy, higher energy 



neutrons are significantly peaked in the forward direction. 
Moreover, the charged particle array has more material at 
forward angles and the Neutron Ball is not symmetric with 
respect to the target. Therefore, the neutron efficiency was 
calculated as a function of both neutron energy and polar 
angle. Further details can be found in Ref. [12]. This 
energy and angle dependent efficiency was implemented in 
the experimental filter. 

The experimental filter and the QP reconstruction pro- 
cedure (i.e. VCut, ZCut and QCut) were applied to the 
simulated data. The reduced impact parameter distri- 
bution of fully equilibrated quasi-projectile fragmentation 
events ranges from 0.5 to 0.9 and it is peaked around 0.7 
for the CoMD calculation. For HIPSE results, the dis- 
tribution is peaked around 0.6 and extends from 0.5 to 
slightly above 0.8. For this reason we will constrain our 
analysis to 0.5 < b re d < 0.9 for the rest of this paper. 
Data were divided into 4 bins, according to the reduced 
impact parameter: b rec ii — 0.55 ±0.05, b re d2 — 0.65 ±0.05, 
b red3 = 0.75 ± 0.05, b red4 = 0.85 ± 0.05. HIPSE-SIMON 
Znresults are not shown for bins which contain fewer than 
100 events. 

3.1. Overall, QP and QT efficiencies determination 

The efficiency of the Neutron Ball as a function of en- 
ergy and angle from the target position for isotropic emis- 
sion is well established [13] ■ However, we are interested 
in reconstructing QPs resulting from multifragmentation 
events. Therefore it is important to know the effective 
neutron efficiencies for these types of events. 

3.1.1. Overall efficiency 

The overall Neutron Ball efficiency was determined as: 



N, 



del 



(2) 



where N raw is the number of neutrons predicted by the 
simulations and Ndet is the number of neutrons passing the 
experimental filter, i.e. the detected neutrons. The calcu- 
lation was performed on both HIPSE-SIMON and CoMD 
simulated data. 

The overall efficiencies, e tot , are shown as horizontal 
lines in Fig. [T] The dotted line corresponds to values ob- 
tained by HIPSE-SIMON calculation for the 70 Zn+ 70 Zn 
reaction, while the other lines correspond to values ob- 
tained by CoMD simulations for the three different reac- 
tions. The calculated uncertainty in the efficiencies re- 
ported throughout the rest of this text are of the order of 
or below 1% and will be omitted in the following. Efficien- 
cies of about 69% are obtained for CoMD simulated data 
for the three reactions. The values obtained for the 64 Zn 
and 64 Ni reactions overlap. No significant dependence on 
the analysed reaction is observed, therefore only values 
obtained for the 70 Zn+ 70 Zn reaction will be shown in the 
remainder of the figures in the paper. A slightly higher e tot 
of 69.2% is obtained from the HIPSE-SIMON simulation 



3 







CoMD 






HIPSE 


0.71 


' — 


• 70 Zn+ 70 Zn 

-■ • 64 7nj- 64 7n 

m £-11+ £-\\ 




□ 


70 Zn + 70 Zn 


0.7 




- • 64 Ni+ 64 Ni 






• 


o 






□ 




□ • 

* 


0.69 














□_ _ _ 


• 










• 
• 

• 


• 
• 






0.68 


T , 






i , 




0.4 


0.5 0.6 




b red 


0.8 0.9 1 



Figure 1: (Color online) Overall efficiencies determined from CoMD 
and HIPSE simulated data. Lines correspond to overall efficiencies 
obtained without sorting the data by reduced impact parameter. 
Lines corresponding to 64 Zn+ 64 Zn and 64 Ni+ 64 Ni overlap. 



(dotted line) for the 70 Zn+ 70 Zn reaction, which is within 
0.3% of that obtained from the CoMD results. 

Overall efficiencies, sorted by reduced impact parame- 
ter, are also plotted for the three analysed reactions in Fig. 
[1] Full and open symbols correspond to values obtained 
from CoMD and HIPSE-SIMON results, respectively. No 
HIPSE simulated data are available in the highest b re d 
bin, due to lack of statistics. Although we observe a slight 
dependence of the total efficiency on the reduced impact 
parameter, the obtained values are in overall agreement to 
within less than 1%. As for the global efficiencies, there 
is no dependence on the reaction and the values obtained 
from CoMD and HIPSE agree to within 1%. 

3.1.2. Efficacy of the velocity selection 

The Neutron Ball efficiency for quasi-projectile emitted 
neutrons was determined as 



N. 



SQPr. 



QPn, det 



QPn, raw 



(3) 



where NQp n ^ raw is the number of QP-neutrons predicted 
by the simulation and NQp 7li det is the number of the de- 
tected QP-emitted neutrons, i.e. those neutrons emitted 
from the QP which pass the experimental filter. QT effi- 
ciencies can be obtained from Eq. [3] by considering QT- 
emittcd neutrons rather than QP. In the HIPSE-SIMON 
output, neutrons can be assigned to the QP based on the 
source tag provided by the model. Also, a velocity con- 
straint, similar to that applied to charged particles, can 
be applied to assign neutrons to the QP. The efficacy of 
such a constraint is investigated by comparing the efficien- 
cies for QP-emitted neutrons obtained using the two selec- 
tions. This is important for the following analysis because, 
in CoMD results, information on the emitting source (i.e. 
the source tag) is not available. 

As discussed in Sec. 12.21 the velocity selection allows 
us to select QP-cmitted charged particles based on their 



longitudinal velocity with respect to the largest fragment 
emitted in the event (projectile-like fragment, PLF). From 
simulation rersults, the selection limits of the neutron lon- 
gitudinal velocity relative to the largest fragment were de- 
termined by correcting for the Coulomb boost that acts on 
the protons. The QP reference frame was reconstructed by 
taking into account all the charged particles accepted by 
the VCut. Neutrons are not considered at this point since, 
experimentally, neutron velocities are not measured and, 
therefore, cannot be included in the QP velocity (vqp) 
determination. The Coulomb boost, Ec, in the quasi- 
projectile reference frame was determined as the difference 
between the proton and neutron average kinetic energies 
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where v$ P ^ and v^ P ^ arc the neutron and proton ve- 
locities in the QP reference frame, respectively. Next, we 
bring Eq. [5] into the lab reference frame and divide by the 
PLF velocity (vplf)- We substitute the used upper and 
lower limits for the quantity Vp ab ^ /vp£ F (1.65 and 0.35, see 
Sec. 12. 2p and assume that the QP velocity is similar to the 
largest fragment velocity (vqp ~ vplf)- This leads to the 
conclusion that the longitudinal neutron velocity relative 
to the largest fragment limits are ±49%. The calculation 
was performed for 64 Zn+ 64 Zn, 70 Zn+ 70 Zn and 64 Ni+ 64 Ni 
reactions, sorting the data by reduced impact parameter. 
The values obtained for the three reactions agree to within 
about 2%. 

It is important to ensure that the the source tag and 
the velocity selection for associating free neutrons with 
the QP produce similar results. For this we investigate 
the percentages of HIPSE-tagged QP neutrons which also 
pass the Vcut. The obtained values are reported for each 
b re d bin in Tab. [T]for the 70 Zn-|- 70 Zn reaction. More than 
90% of QP-tagged neutrons pass the velocity selection for 
bred > 0.6, and about 85% are accepted for the most cen- 
tral considered b re d range. A similar analysis was per- 
formed on QP-tagged protons and the obtained values are 
also reported in the table. A very good agreement is ob- 
served between the values obtained for neutrons and pro- 
tons. In both cases the velocity selection is very effective 
in accepting QP-emitted neutrons and protons, and the se- 
lection only gets better for the more peripheral collisions. 

The capability of the VCut to select only QP-emitted 
neutrons (and protons) is investigated by looking at the 
emitting source of particles accepted by the velocity selec- 
tion (i.e. FCW-acccpted particles). Sources will be clas- 
sified as QP, QT and "other sources", the latter includ- 
ing emitting sources other than the QP and QT. Table [5] 
shows the emitting source, according to HIPSE, of VCut- 
accepted neutrons and protons for the 70 Zn-r- 70 Zn reaction. 
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bred 


0.5-0.6 


0.6-0.7 


0.7-0.8 


0.8-0.9 


neutrons 


85.5% 


93.2% 


96.7% 


98.5% 


protons 


82.6% 


90.3% 


95.4% 


98.6% 



Table 1: Percentages of QP-tagged neutrons and protons accepted 
by the velocity selection for each reduced impact parameter bin for 
the 70 Zn+ 70 Zn reaction, as predicted by HIPSE-SIMON simulation. 
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Figure 2: (Color online) Fraction of VCui-accepted neutrons (full 
symbols) and protons (open symbols) produced by each source rela- 
tive to the total number of neutrons that pass the VCut, according 
to HIPSE-SIMON calculation. Full and dashed lines are the values 
obtained without sorting the data by b re( j, for neutrons and protons 
respectively. 



VCwi-accepted neutrons (and protons) mainly (w 70%) 
come from the QP source. Quasi-target and "other source" 
contributions are of the order of « 10% and 20%, re- 
spectively. These values are plotted in Fig. [2] as full (emit- 
ted neutrons) and dashed (emitted protons) red, blue and 
black lines for QP/QT/"other source". Values sorted by 
b rec i are also plotted in Fig. [5] as circles, squares and tri- 
angles for QP, QT and "other sources". Open and full 
symbols correspond to protons and neutrons, respectively. 
The fraction of VCut- accepted neutrons belonging to the 
QP source (red circles) increases as the b rec i increases while 
the contamination from neutrons emitted by the QT and 
"other sources" decreases, consequently. Similar results 
are obtained for the other reactions. The velocity selec- 
tion appears to be very effective in selecting QP-neutrons 
throughout the range of b re d and gets increasingly better 
for more peripheral collisions. 

Quasi-projectile efficiencies were calculated according 
to Eq. [31 using both QP-tagged and FCtti-accepted neu- 
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QP 
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NV Cut (source)/NV Cut 

N VCut( SOUrCe )/ N VCut 


70% 
68% 


12% 
12% 


18% 
20% 



Table 2: Percentage of VCut-accepted neutrons and protons pro- 
duced by each source relative to the total number of neutrons pass- 
ing the VCut for the 70 Zn+ 70 Zn reaction, as predicted by HIPSE- 
SIMON simulation. 
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75.6% 
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Table 3: Quasi-projectile (egp) and quasi-target (eqt) efficien- 
cies calculated according to Eq. [3] for both models. Tagged and 
VCui-accepted neutrons were used to determine £Qp(tagged) an d 
ZQP(Vcut-acc), respectively (see text). 



trons from HIPSE-SIMON. The obtained e QPn , integrated 
over the 0.5 — 0.9 reduced impact parameter range, are 
reported in Tab. [3] for all reactions as £Qp(tagged) an d 
£ QP(Vcut-acc), respectively. The obtained values are in 
good agreement and no dependence on the reaction is ob- 
served. In Fig. [21 QP neutron efficiencies are plotted for 
the 70 Zn+ 70 Zn reaction as a function of the reduced im- 
pact parameter as determined by QP-tagged neutrons (cir- 
cles) and by velocity-accepted neutrons (squares). The av- 
erage quasi-projectile neutron efficiencies are also plotted 
for QP-tagged neutrons (full line) and velocity-accepted 
neutrons (dashed line). A small dependence on the im- 
pact parameter is observed. Nevertheless, the obtained 
values arc in agreement. Such agreement indicates that 
the velocity selection is an effective way to identify neu- 
trons emitted from the quasi-projectile source. 

From an experimental point of view, the velocity selec- 
tion cannot be applied to neutrons. In order to reconstruct 
a QP, then, the charged particles and neutrons must be 
assigned to the QP differently - by the VCut for charged 
particles and by Eq. [1] for neutrons. In the simulation, 
our newly derived neutron velocity selection criteria be- 
haves very similarly to the proton velocity cut that has 
been previously used to reconstruct QP's and, therefore, 
suggests that it is a successful procedure for assigning neu- 
trons to the QP. Furthermore, this neutron velocity selec- 
tion is consistent with assignment of neutrons to the QP 
based on the source tag from HIPSE. We conclude, then, 
that, within the results of our simulations, the method of 
QP neutron assignment based on the source tag selection 
method, which is unavailable for CoMD results, is an ef- 
fective substitute for the velocity cut selection. 

3.1.3. Quasi-projectile and quasi-target efficiencies 

Quasi-projectile and quasi-target efficiencies were cal- 
culated according to Eq. [3] by considering QP and QT- 
tagged neutrons from the HIPSE simulation and by con- 
sidering FCW-accepted neutrons from CoMD simulation. 
The number of QP-emitted neutrons predicted by the CoMD 
simulation {NQp n ^ raw ) was determined on the basis of the 
longitudinal neutron momentum. Indeed, the distribution 
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Figure 3: Quasi-projectile neutron efficiencies calculated according to 
cq[3]using QP-tagged {e QP ^ agged) ) and KCui-accepted (£Q P (y cut )) 
neutrons from HIPSE-SIMON results for 70 Zn+ 70 Zn reaction. Full 
and dashed lines are the values obtained using all b re d (see text). 
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Figure 5: (Color online) Multiplicity of emitted free neutrons {M 1 n aw ) 
versus multiplicity of detected free neutrons (M^ et ) as obtained from 
HIPSE-SIMON (a) and CoMD (b) simulations for the 70 Zn+ 70 Zn 
reaction (see text). 



Figure 4: (Color online) Quasi-projectile (circles) and quasi-target 
(triangles) emitted neutron detection efficiencies obtained for the 
70 Zn+ 70 Zn reaction from CoMD (full symbols and line) and HIPSE- 
SIMON (open symbols and dashed line) simulated data. 



of neutron momenta presents two peaks that are clearly 
separated, especially for the most peripheral b re d- These 
peaks correspond to neutrons emitted from the QP (for- 
ward in center-of-mass) and the QT (backward in center- 
of-mass) sources. The QP/QT efficiencies are reported 
in Tab. [3] for all the reactions and for both the models. 
Values of about 76% and 61% arc obtained for eqp and 
sqt, respectively No dependence on the system is ob- 
served. Moreover, the values obtained for HIPSE-SIMON 
and CoMD simulated data are in agreement within less 
than 2%. This suggests that these values arc largely model 
independent. 

In Figdl eqp (black lines) and sqt (red lines) are 
plotted for the 70 Zn+ 70 Zn reaction for CoMD (full lines) 
and HIPSE (dashed lines) simulations. The efficiencies for 
neutrons sorted by impact parameter are also plotted in 
FigdJ Full and empty symbols indicate values obtained for 
CoMD and HIPSE-SIMON simulated data, respectively. 



The efficiency value obtained from HIPSE results for the 
most peripheral collisions (b re d = 0.8 — 0.9) is affected by 
a larger uncertainty due to lack of statistics. The CoMD- 
extracted efficiency value for b re d = 0.5 — 0.6 is not pre- 
sented because for the more central collisions the distribu- 
tion of neutron momenta does not present two separated 
peaks. No significant dependence on b re d is observed and 
there is good agreement between the values obtained from 
CoMD and HIPSE-SIMON simulations. 

3.2. An event-by- event analysis of the Neutron Ball 

Once the overall, QP and QT free neutron detection ef- 
ficiencies for multifragmentation events were determined, 
we investigated the event-by-event response of the Neu- 
tron Ball. In Fig. [5] the multiplicity of free neutrons for 
each event (M™ 1 ") is plotted versus the multiplicity of de- 
tected free neutrons (M^ et ), i.e. the number of neutrons 
that pass the experimental filter. Figure O^a) and 0Jb) 
were obtained by analyzing the HIPSE-SIMON and CoMD 
simulated data, respectively, for the 70 Zn-|- 70 Zn reaction. 
The two models predict, for the same reaction, very dif- 
ferent numbers of emitted neutrons. Nevertheless a rather 
narrow correlation of MI at0 and M!^ et is observed in both 
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Figure 6: (Color online) Multiplicity of QP-cmitted free neutrons, 
M^ aw (unfiltered), versus multiplicity of detected QP-emittcd free 
neutrons, M% et (filtered), as obtained from HIPSE-SIMON (a) and 
CoMD (b) simulations for the 70 Zn+ 70 Zn reaction. Neutrons are 
assigned to the QP source according to the tag and the VCut for 
HIPSE-SIMON and CoMD calculations, respectively (see text). 



cases. The standard deviation, a, of the total multiplicity 
M™ w is approximately constant for each M^ et value for 
both models. The weighted average a is 1.92 and it varies 
from 1.5 to 2.0 for < M* et < 17 for HIPSE results. For 
CoMD simulated data the average a is 2.5 and it varies 
from 1.2 to 3 for < Mff < 26. The number of detected 
neutrons, M^ et , is smaller than the number of raw neu- 
trons, consistent with a detection efficiency smaller than 
1. The observed correlation shows that information on 
the real free neutron multiplicity can be extracted from 
the experimentally measured neutron number. 

The effect of the experimental filter on the number of 
QP-emitted neutrons was also investigated. In Fig. [[5] the 
total number of QP-emitted neutrons is plotted versus the 
number of detected QP-emitted neutrons. Figure EJa) is 
obtained for HIPSE simulated data and QP-tagged neu- 
trons are considered. The narrow correlation observed for 
the raw neutron multiplicity and the multiplicity of de- 
tected neutrons (Fig. 0) is observed also between the raw 
QP-emitted neutron multiplicity (Mq^) and the multi- 
plicity of detected neutrons emitted by the QP (Mgp^). 



In this case the width is slightly smaller, with an average 
value of (a) = 1.22, and does not vary significantly with 

M det 

Figure^b) is obtained for CoMD simulated data. Neu- 
trons are assigned to the QP source according to the ve- 
locity selection described in Sec. 13.1.21 Once more, the 
correlation between Mq™ and Mqp^ is observed, with a 
constant average width of 2.19 neutrons. From these ob- 
servations we can conclude that the selection of the neu- 
tron emission source does not wash out the strong corre- 
lation that we see in both panels of Fig. \5\ The Neutron 
Ball response remains sharply peaked which suggests re- 
liable event-by-event neutron multiplicity measurements. 
Moreover, it speaks to the quality of the velocity selection 
that, when applied to CoMD results, the correlation does 
not broaden. Similar narrow correlations are observed for 
QP-emitted protons and for each of the other reactions. 
This gives us confidence that, with a proper treatment of 
the efficiencies, we are able to determine the number of 
neutrons to assign to the QP. 

3.3. Accuracy and precision of the QP neutron assignment 
method 

Knowing the source-specific efficiencies for multifrag- 
mentation events, as determined from HIPSE or CoMD 
simulations, allows one to calculate the number of QP- 
emitted neutrons from the total number of detected free 
neutrons using Eq. [TJ We recall that eqp and eqt were 
found to be model independent in our analysis. The total 
number of detected neutrons (M^ et ) predicted by HIPSE- 
SIMON and CoMD calculations were used, event-by-event, 
to calculate the number of neutrons to be assigned to the 
QP (Nqp). In FigJ7]the raw number of QP-emitted neu- 
trons (Mqp^) is plotted against the multiplicity of QP- 
assigned neutrons as predicted by HIPSE-SIMON (a) and 
CoMD (b) simulations. The full line indicates Nqp = 
Mgp". Each Nqp value corresponds to a distribution of 
raw QP-emitted neutrons predicted by the simulation. In 
order to better depict the correlations in Fig. [7] we can 
obtain a distribution of raw QP-emitted neutrons for each 
bin of calculated QP neutrons. In Tab. U we report the 
mean and widths of those distributions for the 70 Zn+ 70 Zn 
reaction. For the HIPSE-SIMON simulated data, Eq. □ 
gives the expected number of neutrons to associate to the 
QP to within 1 neutron. This is true for events with be- 
tween 3 and 6 neutrons added back into the QP, which cor- 
respond to about 52% of the analysed events. For events 
with Nqp = — 8, the calculated Nqp value is in agree- 
ment with the raw QP-emitted neutron multiplicity within 
less than 1.5 neutrons. This corresponds to about 92% of 
the analysed events. The expected number of neutrons 
to associate to the QP is reproduced by Eq[T] for CoMD 
results to within 1 for about 70% of the analysed events. 
This corresponds to events with 6 or fewer neutrons added 
back to the QP. For events with Nqp = 0—11, which cor- 
respond to about 92% of the available statistics, the calcu- 
lated Nqp is in agreement with the raw QP-emitted neu- 
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tron multiplicity within less than 1.5 neutrons. Therefore, 
this method is reasonably accurate and sufficiently precise 
as to allow for the study of well-defined QP sources. Nev- 
ertheless, the width of the distribution of the raw number 
of QP-emitted neutrons for each value of the calculated 
multiplicity of QP-assigned neutrons introduces an uncer- 
tainty into the determination of the QP composition. In 
any analysis, the impact of such uncertainty on the results 
should be taken into account. 
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Table 4: Average {{Mgp^ )) and a of the QP-emitted neutron multi- 
plicity distribution, plotted in Fig[7] for each value of reconstructed 
QP neutron multiplicity {Nqp) for HIPSE-SIMON and CoMD sim- 
ulated data. 
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Figure 8: (Color online) Reconstructed experimental QP neutron 
multiplicity as given by Eq[T] for the 70 Zn+ 70 Zn, 64 Ni+ 64 Ni and 
64 Zn+ 64 Zn systems. 



Figure 7: (Color online) Multiplicity of QP-emitted neutrons 
{M^Qp ) versus multiplicity of free neutron assigned to the recon- 
structed QP, Nqp, according to Eq[T] for the 70 Zn+™Zn. (a. 
HIPSE-SIMON and b. CoMD calculations) 



4. Experimental QP neutron distributions 

As an example, the described criteria for assigning de- 
tected free neutrons to the QP fragmenting source was 
applied to experimental data on 70 Zn-r- 70 Zn, 64 Ni+ 64 Ni 
and 64 Zn+ 64 Zn measured with the 4?r NIMROD-ISiS array 
described in Scc l2.ll The QP source was reconstructed, 
event-by-event, applying sequentially the VCut, the SumZ 
(the total detected charge was required to be within 25 
and 30), and the QCut to the data. To be applied to 



experimental data, Eq. [T] needs to be modified as: 

N det ~ N b 

ackground 



N, 



QP 



(6) 



The measured neutron multiplicity (Ndet) is corrected for 
the measured background multiplicity {Nb ac kg ro und) > which 
is determined using the gating system discussed in Sec. 

o 

The £cf I^CfoEANT term is the ratio of the efficiencies 
for a 252 Cf source positioned at the target position de- 
termined during the experimental campaign (scf) an< i by 
the GEANT-3 simulation (£c/ G eant)i described in Sec. 
[3] Indeed the GEANT-3 simulation, on which the exper- 
imental filter is based, predicts a smaller number of neu- 
trons to be detected with respect to the experimental one. 
This is due to a not perfect reproduction of all the pro- 
cesses involved in the neutron detection in the simulation. 
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Therefore the efficiencies (eqp and sqt) that GEANT-3 
calculates are lower than the real values. While this is 
not important when dealing with simulated data, since 
both Ndet and the QP and QT efficiencies are reduced of 
the same factor, it has to be taken into account in ex- 
perimental data. This is done by introducing the term 
£ c f /ec fcEANT m Eq. El The experimental scf was mea- 
sured to be 0.7, while the GEANT-3 value was 0.6. As 
for the QP and QT efficiencies, egp and £qt, we used 
the values obtained from the HIPSE-SIMON calculation 
reported in Tab. [3] Equation [5] was applied to the exper- 
imental free neutron multiplicities measured by Neutron 
Ball for the three reactions and allowed the determina- 
tion of the QP mass. The QP N/Z distributions are cen- 
tered around 1.35, 1.30 and 1.20, with a width (RMS) of 
0.13, for the 70 Zn+ 70 Zn, 64 Ni+ 64 Ni and 64 Zn+ 64 Zn reac- 
tions, respectively, as it was shown in In Fig. [5] we 
present the obtained reconstructed QP neutron multiplici- 
ties for the 70 Zn+ 70 Zn (circles), 64 Ni+ 64 Ni (triangles) and 
64 Zn+ 64 Zn (squares). The distributions were normalized 
to 1. The highest mean Nqp value {(Nqp) = 7.95) is ob- 
tained for the most neutron-rich reaction ( 70 Zn+ 70 Zn) and 
it decreases as the average QP isospin decreases. Indeed 
(Nqp) values of 6.69 and 5.52 are found for 64 Ni+ 64 Ni 
and 64 Zn+ 64 Zn, respectively. These observations are con- 
sistent with the expected enhanced emission of neutrons 
for more neutron-rich systems. 

5. Conclusions 

The QP reconstruction procedure - in particular the 
method used to determine the multiplicity of QP-cmitted 
free neutrons - was investigated in detail through two dif- 
ferent simulation codes, CoMD and HIPSE-SIMON. 

The ability of the velocity cut to select QP-emitted 
neutrons (and protons) was carefully investigated. The 
methods based on the velocity selection (typically used 
for charged particles) and on efficiencies calculated from 
QP and QT-tagged neutrons from HIPSE are consistent. 
Both CoMD and HIPSE-SIMON simulated data were used 
to determine the overall and source-specific detection effi- 
ciencies of the TAMU Neutron Ball for multifragmcntation 
events. The agreement within 2% of the results obtained 
from the two codes suggests the efficiency values to be 
largely model-independent. The narrow correlation ob- 
served between the actual emitted and detected neutron 
multiplicities indicates that measured free neutron multi- 
plicities can be used to reproduce the actual free neutron 
multiplicities, even on an event-by-event basis. The sim- 
ilar narrow correlations observed for predicted and mea- 
sured QP-emitted neutrons suggest that, with a proper 
treatment of the efficiencies, the number of neutrons to be 
assigned to the QP can be determined with an accuracy 
of around 1.5 neutrons on average. 

Finally, we showed that, with the efficiencies obtained 
from our calculations, the method using Eq. [T] to exper- 
imentally determine the QP-emitted free neutron multi- 



plicity is reasonably accurate and sufficiently precise as to 
allow for the study of well-defined QP sources. Neverthe- 
less, there remains an uncertainty on the QP composition, 
which should be taken into account by any analysis relying 
on the QP composition. The model independence of our 
results of the analysed reaction suggests that the described 
procedure using Eq. Q] and the obtained efficiency values 
can be used in future analysis with different systems and 
energies provided that efficiencies are recalculated for the 
different reacting systems. 
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